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We investigate the influence of walls and corners (with Dirichlet and Neumann boundary conditions) in 
the evolution of twodimensional autooscillating fields described by the Complex ginzburg-Landau equation. 
Analytical solutions are found, and arguments provided, to show that Dirichlet walls introduce strong selection 
mechanisms for the wave pattern. Corners between walls provide additional synchronization mechanisms and 
associated selection criteria. The numerical results fit well with the theoretical predictions in the parameter 
range studied. 



I. INTRODUCTION. 

Spatially extended nonlinear dynamical systems display an 
amazing variety of behavior including pattern formation, self- 
organization, and spatio-temporal chaos pl-fl]. Much effort 
has been devoted to the characterization of the different dy- 
namical regimes and the transitions between them for model 
equations such as the complex Ginzburg-Landau equation 
(CGLE) [j|]. This is an equation for a complex field j4(x, t) 
that conveniently rescaled reads 



t A = A + (1 + ia)V 2 A - (1 + il3)\A\ 2 A 



(1) 



a and (3 are real parameters. This equation describes the onset 
of an oscillatory regime through the Hopf bifurcation of a ho- 
mogeneous state, and it is used generally as a model equation 
due to the rich variety of its solutions. Binary fluid convection 
, transversally extended lasers [|8j^], chemical turbulence 
jUt], bluff body wakes the motion of bars in the bed 
rivers fll3|], and other systems |Jl4|] have been described using 
the CGLE in a proper parameter range. 

The CGLE admits simple plane- wave solutions. How- 
ever, for most of the (a, (3) parameter range, a typical evolu- 
tion starting from random initial conditions leads to complex, 
steady or evolving, states. An important ingredient in the de- 
scription of these dynamical regimes in two-dimensional do- 
mains is the interaction of singular points called defects. For 
our purposes, a defect is just a zero of the complex field A, 
where there is a singularity in the phase <p defined by the re- 
lation A = \A\ exp(z<z>). There is a topological charge associ- 
ated to each defect, n, defined by 



— </" V(j3 • dr 



(2) 



where V is a closed path around the defect. The topological 
nature of the phase singularities implies that n is a positive 
or negative integer, and that the total topological charge in 
the two-dimensional system is constant, except for the defects 
flowing in and out through the boundaries. In the interior of 
the system, defects can only be created or annihilated by pairs 
of opposite charge. Spiral defects, i.e. defects around which 
the lines of constant phase have a spiral form, are typically 



formed in the CGLE (for ct / The interaction between 
these spiral structures has attracted much attention [0]. Spiral 
solutions of a different nature appear e.g. in excitable media 
such as the Belousov-Zabotinsky reactio n Jlg jlq] and electro- 



hydrodynamic convection (see e.g., ref ||17j|) 

One important source of defects in real systems are the 
boundaries. Under some circumstances, walls can introduce 
defects into the system increasing the amount of disorder in 
the dynamics. In other situations the boundaries play the op- 
posite role: they annihilate defects driving the system to a 
more ordered state. In general, the interplay between these 
two behaviors and the proper dynamics of the bulk can push 
the system to configurations different from the ones found un- 
der boundary-free conditions (periodic boundary conditions 
for instance). However, few studies have been addressed to 
the influence of the boundary shapes and boundary conditions 
on complex dynamics. The importance of these effects in the 
transverse patterns of laser emission, where aspect ratios are 
not large, is visible in recent works such as Jl^ , p^| ]. In ad- 
dition, average patterns in Faraday waves and other spatio- 
temporally chaotic systems have been observed to be sensible 
to boundary shape [ p0pl| ] and boundaries are able even to in- 
duce spatial chaos in otherwise non-chaotic systems [p2j]. All 
those strong influences of boundaries on the dynamics of ex- 
tended nonlinear systems [^3j provides us with the motivation 
for a more systematic study of boundary effects on the CGLE 
performed in this Paper. 



In Ref. [24 1 we performed a first numerical exploration of 
these effects, via computer simulations of the CGLE in cir- 
cular and rectangular geometries with null Dirichlet boundary 
conditions. The results reveal a fundamental role of bound- 
aries in selection mechanisms. In particular wave emission 
from Dirichlet walls (i.e., walls where A = 0), and the dom- 
inance of corners as pacemakers for the whole system were 
important observed effects. Understanding the origin of such 
effects is the main goal of this Paper. To achieve it, we will 
focus first on the effect of a single lateral wall, where the com- 
plex field is set to zero, in the selection of the pattern. After 
this we will study how the presence of corners (i.e. the in- 
tersection of two walls) restricts the family of solutions found 
previously. It should be noted that we use Dirichlet boundary 
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conditions (and in some places of this Paper also Neumann 
boundary conditions) as simple phenomenological conditions 
to explore deviations with respect to the more commonly used 
periodic boundaries. A different issue is to establish what are 
the pertinent boundary conditions arising for the CGLE when 
it is derived as an amplitude equation in particular physical 
contexts (for example in optics, fluids, etc.). Some results in 
this last subject can be found in Ref. [HQ. 

In the next Section we review previous numerical results on 
the CGLE in several geometries and boundary conditions. In 
Section |e[ we summarize analytical solutions in unbounded 
domains. In Section IV, we present analytical and numerical 
results for the CGLE in the presence of a lateral Dirichlet wall. 
In Section M, we extend our study to the case of corners and 
in Section Wu we finish with our Conclusions. 



II. NUMERICAL OBSERVATIONS. 



possibly complex function of position. 

In the final frozen state, defects are confined in the center 
of the domain forming a rigid static chain. Shock lines ap- 
pear where waves from different sides of the contour collide. 
The strongest shocks are attached perpendicularly to the walls 
and the general shock configuration is what one would expect 
for small symmetry breaking of the square geometry p7[ ] The 
number of defects depends on the initial condition. Solutions 
with no defects are also found (e.g Fig. ||(c,d)), and are called 
target-like solutions. This kind of solutions is not seen in sim- 
ulations with periodic boundary conditions. 

In our simulations in the square geometry with Dirichlet 
boundary conditions, the direction of the phase velocity (from 
the walls or towards the walls) and the wavenumber depend 
on the parameter values in a way which differs from the usual 
spiral waves in infinite systems (see Ref. |^] and Section |li|). 
Thus boundaries are playing an important role in the selection 
of the wave speed and wavenumber. 



It is quite evident, and confirmed by our previous study 
p4|], that the effect of boundaries is more noticeable in the 
parameter regimes for which large correlation lengths are 
present in the system. In strongly chaotic states with short 
correlation lengths, the main effects of walls are restricted to 
boundary layers close to them. In consequence we restrict 
here the presentation of our numerical results to the region of 
parameters for which coherent oscillations extend over nearly 
the whole system. This happens in most of the Benjamin-Feir 
stable region in parameter space, that is, for 1 + a/3 > 0, but 
also in other regions close to it. Defects and shocks however 
disrupt the otherwise ordered plane waves, and its location and 
structure are strongly dependent on boundaries. In Fig. |IJ the 
CGLE is solved in a square with null Dirichlet boundary con- 
ditions (A = 0). The zero-amplitude boundaries facilitates 
the formation of defects near the walls. Starting from ran- 
dom initial conditions, defects are actively created in the early 
stages of the evolution (See Fig. [l]). After some time how- 
ever all the points on the boundaries synchronize and oscillate 
in phase so that plane waves are emitted. Defect formation 
ceases, and the waves emitted by the walls push the remain- 
ing defects towards the central region of the domain. There 
the defects annihilate in pairs of opposite charge and as a re- 
sult of this process a bound state is formed by the surviving set 
of equal-charge defects. The orientation of the waves emitted 
by the boundaries also changes during the evolution. The syn- 
chronized emission of the early stages produces wave propa- 
gation perpendicular to the boundary but in the late states the 
wavevector tilts to some emission angle of approximately 45 
degrees. The precise value of this angle depends on both the 
parameter values and the geometry of the boundaries. The fact 
that this angle is not exactly 45 degrees is made evident by the 
slight mismatch between the waves coming from orthogonal 
walls. Finally the system reaches a frozen state of the type 
displayed in Fig. ^. The term frozen is used here to denote 
that the modulus is a steady solution, although the phase is 
time-periodic. More concretely, our frozen configurations are 
well-described by A(x) = f(x)e~ luJt , with li real and / a 




FIG. 1. Time evolution of the solution of the Eq. (|l|) at arbitrary 
times with parameter values a = 2, f3 = —0.2. The initial condi- 
tions is random. When the boundary starts emitting waves, the spiral 
defects are pushed to the interior of the domain and annihilate by 
pairs of opposite charge. The modulus of the field is plotted in the 
left column and the phase in the right, in grey-scale. The final state 
(not shown) contains a single defect, as the one in Fig. En 



In a circular domain (Fig. ||), the frozen structures are ei- 
ther targets (no defects) or a single central defect. Groups of 
defects of the same charge can also form bound states, but 
instead of freezing they rotate together. This contrasts with 
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the behavior of the square domains and is correlated with the 
absence of shock lines linking the boundaries to the center in 
the case of the circular domains. These links are probably re- 
sponsible for providing rigidity to the stationary configuration 
in the square case. Tiny shock lines associated to small de- 
partures from circularity in the lines of constant phase can be 
observed also in the circle but these lines end in the bulk of the 
region before reaching the boundaries. On the other hand, the 
constant-phase lines reach the boundaries nearly tangentially 
in contrast to what we observe in the square. For circular do- 
mains the phase-velocity direction changes with parameters 
in a way more similar to the infinite-system spiral. This is 
another feature revealing that circular boundaries introduce 
less rigidity into the pattern than square ones. The absence 
of corners is probably the main qualitative difference. In fact, 
corners are observed to act as pacemakers from which wave 
emission entrains the whole system [|2-i]|. 




FIG. 2. Frozen structures under null Dirichlet boundary condi- 
tions in a square of size 100 x 100. Parameter values are a = 2, 
P = -0.2 (a,b), and a = 2, j3 = -0.6 (c,d). Snapshots of the mod- 
ulus | A\ of the field are shown in the left column and snapshots of the 
phase in the right column. Grey scale runs from black (minimum) to 
white (maximum). 




FIG. 3. Frozen structures under null Dirichlet boundary con- 
ditions in a circle of diameter 100 for parameter values a = 2, 
(3 = -0.2 (a,b), and a = 2, (3 = -0.6 (c,d). Snapshots of the 
modulus \A\ are shown in the left column and the phase is shown in 
the right column. Gray scale as in Fig. hL 



The stadium shape (Fig. mixes features of the two ge- 
ometries previously studied: it has both straight and circular 
borders. In this case, the curves of constant phase arrange 
themselves to combine the two behaviors described above. On 
the one hand the lines meet the straight portions of the border 
of the stadium with some characteristic angle, as it happens in 
square domains. However, these lines bend to become nearly 
tangent to the semi-circles in the places where they meet with 
these portions of the boundaries. A typical frozen solution 
displays a shock line connecting the centers of the circular 
portions of the domain. This shock line usually contains de- 
fects and their dynamics in this stage is much slower than the 
annihilation that occurs in the bulk of a domain without the 
presence of shocks. It is also possible to find defect-free tar- 
get solutions as in the case of the circle, and the behavior of 
the phase velocity is also similar in the way its direction can 
be changed by modifying the parameters. 
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observations is the goal of the next Sections. 



(a) 



III. SOLUTIONS IN UNBOUNDED DOMAINS. 



(b) 




FIG. 4. Frozen structures under null Dirichlet boundary condi- 
tions in domain with stadium shape of axis 200 x 100 for parameter 
values a — 2, j3 = —0.2. Snapshots of the modulus \A\ are shown 
in the left column and the phase is shown in the right column. Gray 
scale as in Fig. |l| 



To summarize, Dirichlet boundary conditions play a double 
role. On one hand, the walls naturally behave as sources (or 
sinks) of defects. On the other hand, a wall with null Dirich- 
let conditions shows a tendency to emit plane waves that will 
coherently fill the whole system. The interplay between these 
two properties of the boundaries gives rise to interesting be- 
havior. In the case of frozen states, the character of the walls 
as wave emitters dominates. The intersection of two walls (a 
corner) is observed also to emit waves, and the whole sys- 
tem becomes synchronized to this corner emission. In circu- 
lar domains, on the other hand, there are no corners and wave 
selection is definitively dominated by the internal spirals. An- 
other aspect of the dynamical dominance of the walls in the 
square case is that defects form a chain which is anchored 
to the boundaries by a set of shock lines; in a circle, on the 
contrary, the asymptotic state is usually a bound state discon- 
nected from the boundaries. 

Gaining some understanding of aspects of our numerical 



In this Section we review some of the analytical solutions 
of the CGLE in unbounded domains. First we start with plane 
waves, continue with one-dimensional holes and finish with 
two-dimensional spirals. 

The CGLE possesses, among many other solutions, a fam- 
ily of plane wave solutions and solutions containing phase- 
singular points. The plane-wave continuous family is parame- 
terized by the corresponding wave number k. The form of the 
solutions is A = i?exp[i(k • x — cot)], where R = vl — k 2 , 
w(k) = ft - k 2 (/3 - a), and k = |k|. The limit of stabil- 
ity of plane waves is known as the Benjamin-Feir line and 
is given by the curve 1 + a/3 = 0; if this quantity is posi- 
tive, some stable plane wave exists [|28|]; if 1 + a/3 < all 
plane waves are unstable. The limit is given by the stabil- 
ity of the plane wave with k = 0. Stability analysis gives that 
plane waves possessing wave number k in the range [—k e , k c ], 

where k c = \+ap*+ 2(P> are sta ble- The instability is with 
respect to long-wavelength disturbances whose wave vectors 
are parallel to k (Eckhaus instability) [|28[]. It will be useful 
for the future discussion to have an expression for the phase 
velocity of the waves, and of the group velocity of small per- 
turbations on such waves, v p ^ and v gr respectively 



u>(k) t 
v f> = -IT k 



»gr 



-2k((3- a)k 



(3) 
(4) 



k is the unit vector in the direction of k. The expression for the 
group velocity | ^0| ] turns out to be equivalent to the linearly- 
looking expression v gr = W\ tL ui(k), even though u(k) is the 
dispersion relation of nonlinear waves. 

In addition to simple waves, the one-dimensional CGLE 
possess a one-parameter family of solutions for which the am- 
plitude displays a region of local depression. Their analytic 
form was determined by Nozaki and Bekki [CT], and they are 
therefore also referred to as Nozaki-Bekki solutions or holes. 
One member the family is characterized by the value of A 
being zero at a point, called the core of the hole, and asymp- 
totically behaving, at both sides of the core, as a plane wave of 
wave number k. It is worth noting this one-dimensional hole 
solution was also obtained by Hagan [ p9| ] as a sub-product of 
his calculations for two-dimensional spirals. At variance with 
the other members of the Nozaki-Bekki family, this hole does 
does not travel into the system and thus it will be denoted as 
the standing hole. Its analytical expression (choosing the ori- 
gin of coordinates at the hole core) can be written as 



Wh(x, t) = \/\ — k 2 tanh(pa;) exp[i(ip(x) — tot)] , (5) 

where tp is a function of x satisfying 

dip/dx — fctanh(px) , (6) 

(i.e. if> = ipo + (k/p) logcosh(pa;), with i/jq an arbitrary ref- 
erence phase) and u>, k, and p are related according to 
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oj = p-k 2 {f3-a) , 
Spa ' 

{4(/3 - a) + 18a(l + a 2 )}/ 

-{4(/3 - a) + 9 a (1 + a/3)}p 2 + /3 - a 

If a — we can easily get p = 1/ \/2 and 

w = /3(l-fe 2 ), 
p=l/y/2, 

(3 = 



3k 



V2(l- fc 2 ) 



(7) 
(8) 

0. (9) 

(10) 
(11) 
(12) 



Thus /3 and fc have opposite sign ((3k < 0), when a = 0. For 
any value of a and (3, the existence of a defect-like solution 
fixes the value of the asymptotic wave number k and accord- 
ingly that of uj. We mention here that for configurations of 
the frozen type, the solutions with a arbitrary can be obtained 
from the ones with a = by means of a change of variables. 
This fact, which frequently simplifies analysis, is detailed in 
the Appendix. 

The phase and group velocity far from the core for the one- 
dimensional standing hole with a = are 



1(1 -fc 2 ) 3 
v gr = -2k(3 > . 



(13) 
(14) 



Thus the propagation of the phase is towards the core of the 
defect independently of the value of (3. However, the group 
velocity is directed outwards from the core of the defect. Thus 
small perturbations to this solution are expelled away from the 
core. The case of arbitrary a can also be solved numerically. 
Given the parameters (a, (3), the line where the phase velocity 
is zero can be found and it is plotted in Fig. |J The group 
velocity turns out to be always positive (i.e. outwards from 
the core) for the standing hole solutions. 

The two-dimensional spiral wave solutions of the CGLE are 
solutions winding around a defect core (i.e. a phase singular- 
ity). In polar coordinates (r, 9) around the core, they have the 
analytical form [ 29 1 : 



D(r, 9, t) = R(r) exp(i(0 + S(r) - tut)) 



(15) 



This solution represents a phase pattern rotating steadily 
around r = with frequency u> (and frozen modulus). The 
amplitude R is a monotonically increasing function of r, pro- 
portional to r near the origin, and asymptotically approaching 
some value R^, < 1 for large r. The function S behaves 
smoothly in the neighborhood of the origin, taking the form 
S ~ So + Sir 2 . Far from the origin S becomes propor- 
tional to r, behaving as S ~ kr. In this way, in the distant 
region, the isophase lines approach the form of Archimedian 
spirals, converging to plane waves with wave number fc. Thus 
^oo = VI — fc 2 and oj = (3 ~ k 2 ((3 — a). The charge of solu- 
tions of the form Eq. ( p"5| ) is, according to eq. (Q), equal to +1. 
There exists also a negatively charged spiral, with the form 



Eq. (|l5]) but with replaced by —9. In spiral waves, wave 
motion is induced in such a manner as to cause the global 
synchronization of the medium by the defect. 

It is important to notice that, in both one (the standing hole) 
and two dimensions (the spiral solution), imposing the re- 
quirements of zero field at the core, and plane wave behavior 
far from the core, the value of fc gets fixed. Thus fixing the 
parameter values (a, f3), the spiral structure (and the standing 
hole) is unique (except for an arbitrary change in the loca- 
tion of the core). The precise way in which wavenumber, fre- 
quency, phase or group velocities depend on parameter values 
(a, f3) can be found for example in Ref. [|2^|. 



ea 




FIG. 5. Parameter space of the CGLE. Different regions are sep- 
arated by solid lines: BF unstable regime where there are no stable 
plane wave solutions in the infinite system; regions where the phase 
velocity v p h is positive or negative are also shown, and separated by 
additional solid lines for the case in which a single Dirichlet wall 
is present in the system. Star line corresponds to zero phase veloc- 
ity for emission from a corner between two Dirichlet walls spanning 
an angle of 135 degrees; diamond line corresponds to a 90-degrees 
corner. 



IV. SOLUTIONS WITH A SINGLE WALL. 



As a first step to understanding the solutions of the CGLE 
in bounded domains, we study in this Section solutions in the 
presence of a single wall where the value of the complex field 
A is set to zero. We observe numerically that, starting from 
random initial conditions in a bounded domain with a sin- 
gle Dirichlet wall, frozen solutions are reached asymptotically 
(see Fig. ||). In our numerical implementation, the Dirich- 
let wall ( A = 0) is the left one, Neumann boundary condi- 
tions (zero normal derivative of A) are applied to the right 
wall, and the upper and lower limits of the domain are iden- 
tified via periodic boundary conditions. Initially some (spi- 
ral) defects are formed. Typically, the Dirichlet wall starts to 
emit plane waves that push the defects towards the opposite 
boundary until they are all expelled or annihilated. The sta- 
tionary solution is the two-dimensional extension of the one- 
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dimensional standing hole solution described in Section |lTl| 
(that is a continuous line of holes with their cores on the wall: 
W H (x,y,t)=W H (x,t)). 

We can investigate the possibility of more complex solu- 
tions in which the amplitude is independent of the y-direction 
and takes the form of a hole solution in one dimension, but 
with a phase that depends explicitly on both coordinates. We 
study first the case of a = to come back later to the general 
case. 

We look for solutions of the form 



Ww(x, y, t) = yl — k 2 t&nh(px) exp[i(?p(x, y) — tut)] 

(16) 

with uj = (3(1 — k 2 ) and k 2 = k 2 + k 2 . Assuming the form 
tp(x, y) = ip(x) + ip(y), we get the relations 

d x ip(x,y) = k x tanh(pcc) 
d y ip(x,y) = k y 

2p 2 = l- k 2 y 
3k xP = 0(1 - k 2 ) 

and substitution of Eq. ( |l9| ) in Eq. ( |20| ) gives 



(17) 
(18) 
(19) 
(20) 



ik x W 1 



L.2 



V2 



-0(1 ~k 2 ) 



(21) 



Note that if k y = we recover the expression for the one- 
dimensional standing hole solution (in particular we recover 
Eq. (0J)). 

We can perform a similar calculation for the general case 
of parameters a and 0. For a solution of the form Eq. (|T^), 
Eqs. |/7|) and ^ remain valid, and u>, k, and p are related 
according to 



k x — 



0-k 2 (0-a) , 
2p 2 + k 2 -l 



Spa 
-3pk x + a(2p 2 



k 2 x )-0(l-k 2 ) 



(22) 
(23) 
(24) 



In contrast with the selection mechanism for the standing 
hole or the spiral solutions, the presence of the wall does 
not select a unique wave vector but a one-parameter family 
of solutions parametrized by either k x or k y arises instead 
from Eqs. ([17])-(|T]) or (p2|)-(p4l) for given values of a and 
0. Different solutions in the family differ in the direction 
and magnitude of the wavevector k. Different wavevectors 
change the angle of intersection between the lines of constant 
phase and the wall. Figures |^(e-f) are the final state in a nu- 
merical simulation in which the initial condition was close to 
( |T^ ) with k oblique with respect the wall. The displayed state 
is identical (far enough from the Neumann wall) to Eq. ( |T6| ) 
with Eqs. (p2[)-(p4)) thus numerically proving the stability of 
this solution. Different orientations of k can be tested in the 
same way. However, if starting with random initial condi- 
tions we typically find solutions corresponding to the case 



k y = (Fig. ||c-d) that is the simplest two-dimensional ex- 
tension (k y = 0) of the standing hole. The prevalence of the 
k y = solution when starting from random initial conditions 
(of small amplitude) can be understood by analyzing the lin- 
ear stability of the state A = with Dirichlet conditions in a 
single wall limiting a semi-infinite domain. The linear eigen- 
modes are of the form a(x) exp(k y y), and it is easy to see that 
the fastest growing one has k y = so that it will overcome 
the other ones at long times before nonlinear saturation. 



We note that, although solutions ( |l6| ) represent emission at 
an angle with the wall, the analytic expression predicts a small 
boundary layer (of size p _1 ) in which the wavenumber leaves 
its asymptotic orientation to become parallel to the wall, so 
that isophase lines arrive perpendicular to the boundary. This 
is observed in the numerical solutions (see for example Fig. ^f, 
and also the rectilinear walls of Figs. [j], |2|, and |j) thus nicely 
confirming the relevance of the analytical solutions to the ob- 
served configurations. The analytical expressions describing 
wave emission with phase lines parallel or at small angle with 
the walls also provide approximate descriptions of the wave 
emission in the circular geometry. In the square, however, 
conflict between the orientation emitted by neighboring walls 
occurs, and the exact expression ( |l6| ) is appropriate only near 
each wall. The conflict between neighboring walls is resolved 
at long times by emission from the corner, as will be seen in 
the next Section. 



Another important kind of solutions with a single wall is 
the one that appear with Neumann boundary conditions. The 
solutions observed close to the right wall in Fig. ^c-f) are of 
this type. These solutions have been already analyzed in the 
literature, specially in the context of interactions between spi- 
rals. The reason is that a Neumann wall acts as a reflecting 
boundary, so that having a wave impinging into the boundary 
is equivalent to the interaction between two sources of waves 
located symmetrically with respect to the wall [|J|. Despite the 
interest of such solutions, no exact analytical expression has 
been found for them. Analytical solutions have been obtained 
however J32] ] by solving the phase equation that approximates 
the phase dynamics for slow amplitude variations. In agree- 
ment with the numerical observations, the solution presents 
a maximum modulus at the wall (a shock) and the isophase 
lines, straight in the far field, deform when entering a bound- 
ary layer close to the wall to arrive parallel (for normal inci- 
dence) or perpendicular (for tilted far-field incidence) to the 
wall. We will see in the following that these 'tilted Dirichlet 
waves' are of relevance when corners are present. 



It is clear on physical grounds, and confirmed by the analyt- 
ical expressions from the phase approximation, that the wall 
can act as a sink of waves of arbitrary far-field orientation and 
wavelength (the maximum modulus at the shock will adapt 
accordingly). Neumann waves constitute, then, a biparamet- 
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ric family for fixed (a, (3) 

(a) 



(c) 



(e) 




FIG. 6. Modulus (left column) and phase (right column) of the 
solution of Eq. (|j]) for a = 2 and f3 = —0.2 with Dirichlet bound- 
ary conditions for the left boundary, Neumann for the right one, and 
periodic in the horizontal ones. (a,b): Early-time state starting from 
random initial conditions of small amplitude, (c-d): The final asymp- 
totic state. The lines of constant phase travel to the right. Notice 
that although there is a developed spiral defect it disappears through 
the Neumann boundary at long times. (e,f): The long-time asymp- 
totic state from an initial condition consisting of distorted plane wave 
with wavevector oblique to the boundaries. A solution of the form 
Eq. dig), with wavevector close to the initial one, is finally reached. 



broken line with a corner of a definite angle <f>. The bound- 
ary conditions are the following: for the right, upper, and bot- 
tom walls, Neumann boundary conditions (null normal deriva- 
tive). For the left boundary (where the corner is present), null 
Dirichlet conditions. This left boundary is a line that is broken 
forming a variable angle <j>. If this angle is 180°, there exist 
the two-dimensional extension of the standing hole described 
in the previous Section. As the angle decreases the wave is not 
longer plane, and the phase velocity adapts to the new geom- 
etry. The wave fronts may become just slightly distorted from 
straight lines (as in Fig. ^ or display a kink (similar to the 
situation in Fig. ||) depending on a, (3, and (j>. In any case the 
kink is never too strong and departures from straight wave- 
fronts never large. Changing parameters the phase velocity 
may vanish. The locus in parameter space where this happens 
is a 2d surface in the (a, (3, 4>) space. Projections in the (a, (3) 
plane for (j> = 180° (obtained from Eq. (§)) and <\> = 135° 
(from numerical simulation) are plotted in Fig. |[ 6 = 90° 
corresponds to a square and is also plotted in Fig. ||. We do 
not see differences between squares with two or four Dirichlet 
walls. 
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V. SOLUTIONS IN PRESENCE OF CORNERS. 

We now pay attention to the effects induced by the presence 
of corners, i.e. how the solutions ad apt t o the emission of 
two semi-infinite lines. In subsection V A we will show that 



the phase velocity not only depends on the parameters of the 
CGLE a and (3, but also on the angle tfi between the walls 
of the boundary. In subsection V B, we will present solutions 



of the phase equation representing wave collision; they are 
usually called V-solutions. These solutions should be matched 
with the boundaries, which provide selection mechanisms for 
the wave pattern. 



A. Phase velocity dependence on the angle at the boundary. 

We have performed numerical simulations of the CGLE in 
the domain depicted in Fig. ^, where one of the walls is a 



FIG. 7. Domain and boundary conditions. In A ad B, null Dirich- 
let boundary conditions; C, D and E, null Neumann boundary condi- 
tions 



Summarizing, for frozen structures, the presence of Dirich- 
let walls establish a selection mechanism different from the 
associated to the presence of a spiral core in an infinite do- 
main (Section ^l|). When the Dirichlet wall is broken, it is 
seen in the earliest stages of wave-pattern development that 
emission with isophases parallel to the walls is initiated, but 
collision between the waves from different walls arises and 
a distinct final state, with wavenumber, phase velocity, etc. 
fixed by (a, (3, (f) is reached. We now investigate how this 
may happen. 
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B. V-solutions of the phase equation and pattern selection. 

For unbounded domains and for small amplitude modula- 
tions, a phase description of the complex field A can be made. 
The approximate equation for the phase is 



ip = oj + b Q V 2 ip + 61 (Vv?) 2 



(25) 



where 60 = 1 + aft, b\ = {(3 — a), and ojq = —(3. 

We look for solutions <p — (f(x, y, t) representing phase 
waves with non-straight isophase lines. This is what is ob- 
served when Dirichlet waves from different waves interact 
(Figs. Analytic expressions of this type are known for 
the phase equation: the V-solutions [^5j]. We impose differ- 
ent but symmetric wavevectors far from a shock occurring at 
x = y, that is V</? — ► (ki, k 2 ) if x <C y, and V</? — > (k 2 , fci) 
si y ^S> x, thus getting the family of solutions: 

tp(x,y,t) = (uo + h(kf + fcf )) t+ kl + k2 {x + y) 



T- log 

h 



, , h ki - k 2 . . 
cosh I — (a; - y) 



(26) 



The spatial dependence of this solution can be separated in 
terms of the variables u — x + y and v = x — y, and thus 
the phase equation is also separable in u and v. The change 
from (x, y) — > (u, v) is a rotation bringing the shock line 
to one of the axes. After inspection of the derivatives nor- 
mal to the shock, we see that half of these V-solutions can 
be interpreted as tilted waves approaching a Neumann wall at 
the shock position, being the other half just a specular image. 
As for Neumann tilted waves, we have a biparametric family, 
parametrized by ki and k 2 . 

As solutions of the phase equation, the V-solutions are 
strictly valid only far from the boundaries, where the mod- 
ulus of the field remains nearly constant. Matching to solu- 
tions of the form of Eq. ( p"6| ) should be performed close to 
Dirichlet boundaries. We know (Sect. IV ) that for this type 
of boundaries, the two components of the far-field wavevec- 
tor are not independent (Eqs. (^)-(E^). This introduces a 
relationship between ki and k 2 in (|26[). If the shock line 
x = y bisects the angle <fr between two Dirichlet walls, no ad- 
ditional constraints appear from matching to the other bound- 
ary. Thus, one of the parameters in the V-solution, which 
can be taken as the angle between the two waves, is still un- 
determined. From the numerical simulations, it appears that 
this angle becomes determined when the medium is synchro- 
nized by the waves coming from the corner between the two 
walls. We do not have a rigorous argument to demonstrate 
that this is the case, but the following heuristic argument is a 
step towards such a demonstration: Close to the walls a phase 
description is not longer valid, and the modulus approaches 
zero. The solution is of the frozen type, which we write as 
A(x,y,t) = R(x,y) exp(i(ip(x,y) — cot)) with real R, ip, 
and u>. Since this solution should become ( p^ ) far from the 
walls, we immediately find 



Sufficiently close to the corner at x = y = 0, the modulus R 
is small. Writing for it and for ip a Taylor expansion, impos- 
ing symmetry across the x — y line, and substituting into the 
CGLE ([[]), we easily find at the lowest order in distance to the 
corner the following behavior: 



R[x,y) « Bxy 
i/j(x, y) w ip - 



w / 2 
— (x 



y 2 ) 



(28) 
(29) 



Close to the walls the local wavevector is q = Vip = 
— y), so that in the diagonal = — x) with 
modulus qd = 77=^. Far from the corner this wavenumber 
should match the one obtained from the V-solution = 
(ki + ki)/y2. An approximate way of doing this is imposing 
that both wavenumbers become equal at some distance iwa 
from the corner: 
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ki + k 2 

~7T~ 



(30) 



a is an unknown constant of the order of the boundary layer 
size (p -1 ). For given parameter values a and (3 this expression 
gives an extra relationship between ki and k 2 or, equivalent! y, 
between the modulus k and the angle of emission from the 
walls. This, and Eqs. (p^-^4|), completely fixes the solution in 
the presence of a corner. 

Of course, precise numerical values can not be obtained 
since a is unknown, but the previous heuristic argument was 
intended only to illustrate how the presence of the corner re- 
solves the conflict between the neighboring waves, and fixes 
the wave pattern as numerically observed. For situations such 
as the ones depicted in Fig. |]for which the wave fronts remain 
relatively straight, we have k\ ~ k 2 , which can be used as a 
substitute of ( |30| ) to fix the pattern. In fact this is never a bad 
approximation. For example, from a 90° corner, straight and 
symmetric wavefronts indicates wave emission at 45° from 
each wall. We have checked that this is in fact equivalent to 
Eq. ( f30| ) with a 2 — 8. Since this value of a is within the 
boundary layer range, both approaches ((|o|) and fci ~ k 2 ) are 
mutually consistent as can be thought as two different approx- 
imations to the same fact that the corner fixes the wavenumber. 
Assuming k\ = k 2 , we have plotted in Fig. ^| a comparison 
between the results from the numerical simulations and the 
analytical predictions. The agreement is good and confirms 
the relevance of the walls and corners into the wave selection 



lu = lo + b 1 (k 2 + fcf) = -j3 + (/3- a) {kl+kl) (27) process. 
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FIG. 8. Phase of the solution of Eq. (|l|) in grey-scale, 
for parameter values a = 2, ft = —0.2. In (a) the angle 

<t> = tt/2 + tan _1 (l/5); (b) <f> = 3tt/4. 



VI. CONCLUSIONS. 

In this Paper we have presented numerical results on the in- 
fluence of boundaries in the wave pattern selection of a selfos- 
cillatory medium, in parameter regimes in which frozen struc- 
tures are reached. Analytical solutions in the presence of walls 
and corners have been presented and shown to be relevant to 
the numerically observed configurations. The dominance of 
Dirichlet walls, the relative passive rol of Neumann walls, and 
the synchronization properties of corners, are possibly generic 
features which should be found in other selfoscillatory model 
systems. Extrapolation to real experimental oscillatory media 
should be made with care, however, since determining the cor- 
rect boundary conditions applying to the amplitude equation 
associated to a particular medium is a subtle task [E5pq]. 
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FIG. 9. Modulus of k versus parameter /3 (for a square with 
Dirichlet walls, and a = 0) obtained from our theoretical arguments 
(Eqs. (^)-(^) and k x — k y , solid line) and direct numerical simu- 
lation (diamonds). 



Note: In the following address we have made available a 
web page containing simula- 

tions of the CGLE in different geometries related to this paper: 
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APPENDIX A 

The equation resulting from restricting the CGLE Eq. ([!]) to 
frozen solutions of the form A(x) = /(x)e~ iwt , with lu a real 
frequency and / a possibly complex function of the position, 
admits a change of variables [ p9[ ] that transforms the case with 
parameters (a, (3) into the case with parameters (0, (3 1 ). The 
transformation is 







(3- a 
l + afi 
uj — a 
1 + uja 



x = x 



/' = / 



H 


- auj 




fa 2 


1 


-a/3 




- uja 



(31) 
(32) 

(33) 
(34) 



Obtaining a frozen solution (i.e., a function /'(x') and an as- 
sociated frequency a/) at parameters a = and j3' thus allows 
finding corresponding solutions (/, oj) at arbitrary values of a 
and the corresponding f3 = 1 j3 _ + " . This useful relationship 



has been used along this Paper to generate solutions at arbi- 
trary parameters from easier solutions at a = 0. Note that 
if /' contains a factor of the form e ik x , then / will have a 
factor of the form e lk x , with k = k'/\/l — auo' . 



URL: nttp : / /www . nbi . dk/CATS| 



URL: xttp : //www . imedea . uib . es /Nonlinear 
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